The connection between noise and quantum correlations in a double quantum dot 



F. Bodoky 

Kavli Institute of NanoScience, Delft University of Technology, Lorentzweg 1, 2628 CJ Delft, The Netherlands 

W. Belzig 

Fachbereich Physik, Universitat Konstanz, 78457 Konstanz, Germany 

C. Bruder 

Departement of Physics, University of Basel, Klingelbergstr. 82, 4056 Basel, Switzerland 

(Dated: February 2, 2008) 

We investigate the current and noise characteristics of a double quantum dot system. The strong 
correlations induced by the Coulomb interaction and the Pauli principle create entangled two- 
electron states and lead to signatures in the transport properties. We show that the interaction 
parameter <j>, which measures the admixture of the double-occupancy contribution to the singlet 
state and thus the degree of entanglement, can be directly accessed through the Fano factor of 
super-Poissonian shot noise. 



PACS numbers: 72.70.+m,73.23.-b 

I. INTRODUCTION 

Quantum dots, their current and noise properties 1 ^ 
are widely investigated topics in today's nanophysics 
with a number of possible applications, the most revo- 
lutionary being a spin-based quantum computer 3 . Here, 
the spin of an electron in a quantum dot is used as a 
qubit, which serves as the basic building block of a quan- 
tum computer. A necessary ingredient for quantum pro- 
cessing is the possibility to couple at least two qubits 
and to understand and manipulate the various correla- 
tions. This is one of the motivations to study double 
quantum dot o i 5 i6 1 7 1 8 1 9 1 io l n,i2 l i3 _ j n p ar ticular, it is re- 
quired to create entangled electron states by the interac- 
tion of an electron inside the double dot with an electron 
tunneling onto the double dot. Measuring this entangle- 
ment is an important experimental task and theoretical 
suggestions how to probe these states are needed. 

In this work, we will discuss the relation between zero- 
frequency current noise, Coulomb correlations, and en- 
tanglement for the example of a double quantum dot. 
We will make use of the recent observatio n 14 ' 15 : 16 that 
noise measurements in the nonlinear Coulomb-blockade 
regime can be used to obtain spectroscopic information 
of excited levels in a quantum dot. This scheme will be 
adapted to a double quantum dot and we will show that 
the important interaction parameter <p (containing infor- 
mation about the entanglement) 7 -^ enters the Fano factor 
(i.e., the ratio of noise power to current) in this regime. 
This kind of noise spectroscopy can be used to directly 
extract (f> from experimentally observable quantities. We 
extend a previous discussion 17 of the low-frequency shot 
noise of a double-dot system by an analysis of the full 
counting statistics of the transferred charge. This al- 
lows us to identify directly the interaction parameter by 
a noise measurement. One main result is obtained in 
the transport regime, in which predominantly only one 
electron occupies the lowest double dot level due to the 



Coulomb interaction. Thermal activation leads to tun- 
neling events through excited two-electron states which 
result (under the conditions discussed below) in super- 
Poissonian shot noise characterized by a Fano factor 

F(0) = l + </> 2 . (1) 

Hence, a measurement of the noise in this regime allows 
a direct determination of the interaction parameter (j). 

II. MODEL AND THEORETICAL METHODS 

We now turn to a description of the model and 
the method. With maximally two electrons occupying 
the double dot and no magnetic field lifting the spin- 
degeneracy, the following Hund-Mulliken eigenstates of 
the double quantum dot are possible (for a more detailed 
description see 4 ^): the zero-electron state |0), two two- 
fold degenerate one-electron states |+) (symmetric) and 
|— ) (anti-symmetric), four two-electron states with one 
electron per dot, the singlet \S) and the three- fold degen- 
erate triplet \T). There are two other singlet states that 
have a significantly higher energy due to the large on-site 
Coulomb interaction which do not have a physical effect 
in the regime studied (but are included in the numerical 
calculations). 

We consider a (longitudinal) double dot coupled to 
leads, which are modeled as a Fermi sea. The tunnel 
coupling leads to a tunneling amplitude t a , which mea- 
sures the overlap between the orbital state and the lead 
wave function in terminal a. We assume that V is the 
applied bias voltage and that fir, = eV/2 (/j,r = —eV/2) 
are the chemical potentials of the left and right leads. 
There is an additional lead with capacitive coupling to 
the double dot characterized by the gate voltage V g that 
shifts the potential on the dot. When one electron is on 
the dot, it can be in two possible states: its wave func- 
tion is spread either symmetrically or anti-symmetrically 
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over the two dots. The symmetric state has a lower 
eigenenergy, and the energy difference between the two 
one-electron states depends on the interdot tunneling am- 
plitude to, describing the potential barrier between the 
two dots. Including the energy contribution of the gate 
voltage, the eigenenergies of these states can be calcu- 
lated as E± = eV g =F to. Adding a second electron, we 
have first to overcome the gate voltage again, and addi- 
tionally the repulsive Coulomb interaction between the 
two electrons leads to a spin-dependent splitting of the 
eigenenergies. For the singlet state \S) we obtain the 
energy E$ = 2eV g + u\2 — J, and for the (three-fold de- 
generate) triplet state \T) the energy Et — 2eV g + u\2- 
Here, M12 is the (inter-dot) electron-electron repulsion 
and J = At\juH is the Heisenberg exchange parameter 
that characterizes the Heisenberg interaction between the 
two spins, H sp i„ = JSi • S2; ujj is the on-site Coulomb 
repulsion. Because J is positive, the singlet state has a 
lower energy than the triplet state. 

The singlet state can be expressed in terms of the 
creation operators S n , n = ±, which create an (anti- 
symmetric electron state 

|5> = TTT^ ~ ^-^-1) |0) ■ (2) 

The interaction parameter <p describes the competition 
between the kinetic-energy gain and the Coulomb repul- 
sion. In the Hund-Mulliken model <j> = \A+ I^o/^h — 
Ato/un is determined by the inter-dot tunneling ampli- 
tude to and the on-site Coulomb repulsion ur~-. In gen- 
eral, the dependence of <fi on the microscopic parameters 
may have a different form, however we will always use this 
definition for the quantitative plots below. It is impor- 
tant to stress that the results obtained below are largely 
independent on the precise dependence of <f> on the de- 
tails of the quantum dot. Hence, the noise features we 
find can be used as an additional test of the applicabil- 
ity of the model calculation of Ref.— to realistic quantum 
dots. 

We describe transport through the double dot in 
the sequential-tunneling approximation using the Mas- 
ter equation 
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Here, i and j label the available states of the quantum 
dot, including the number of electrons on the dot, the 
orbital index, and the spin. P t is the occupation proba- 
bility of state i, and I\j is the tunneling rate from state 
i to state j. The rates are given by 

Tji = 'YLmijf(a ij (eij - /i L )) + 7j i m !3 /((J t j(e 1 j - hr)), 

Here, jl/r — v L/R\tL/R\ 2 is the bare tunneling rate 
to the left (right) terminal, and v L / R is the density of 
states in the leads. The matrix elements my depend on 
the orbital and the spin state of the dot£. The occupation 



in the leads is determined by the Fermi functions /(e) = 
1 / [exp(e / k^T) + 1], the direction of tunneling is defined 
by aij = ±1 and the energy difference by = E,- L — Ej. 
We note here, that the interaction parameter </> enters 
the calculations through the matrix elements my of the 
singlet state. For a more detailed discussion of these 
tunneling rates we refer tc 8 . For later use we express 
the Master equation ([3J in matrix form, dV/dt = MP, 
where M i3 = T jt for i ^ j and M« = - r y- 

From the Master equation, we obtain the 
nois o 18 ' 19 i 20 ' 21 i 22 . The quantity characterizing the 
current noise is the Fano factor F = S/2eI. It measures, 
how far the distribution of the tunneling events differ 
from the random, Poissonian distribution: for Poissonian 
noise, when the tunneling events are uncorrelated, F = 1. 
A Fano factor F < 1 can be due to anti-correlations (e.g. 
the Pauli principle). We will be mostly interested below 
in super-Poissonian noise, F > 1, and we will show that 
here it is due to correlated transport cycles resulting 
from a blockade of open channels. By a transport cycle 
we mean the following process^: suppose the system 
is initially in its ground state A, and another state B 
(which would be available for tunneling according to its 
energy) is blocked due to Coulomb interaction. When 
this electron in A tunnels out of the dot due to a thermal 
excitation, both the states A and B are available for 
electrons to tunnel into. A cycle is occurring, when a 
sequence of electrons tunnels through state B before 
an electron enters the ground state A and blocks the 
transport through the dot again. This cycle leads to a 
correlated transfer of a number of charges given by the 
number electrons tunneling through the excited state B. 
Hence, that effective charge transferred in one cycle is 
larger than one and leads to an increased Fano factor. 

The above picture is based on the sequential-tunneling 
approximation. The parameters have to be chosen such 
that cotunneling processes can be neglected. This is ex- 
perimentally possible, since the current due to cotunnel- 
ing processes is generally of the order of t 4 / A, where t is 
the tunneling amplitude and A is the excitation energy in 
the virtual intermediate state. In our case the current in 
the most interesting regime is proportional to t 2 x, where 
x is exponentially small in A/k^T. Since the tunneling 
amplitude t can in principle be made arbitrarily small, 
sequential tunneling becomes the dominant process. Re- 
cently it was shown that noise spectroscopy deep inside 
the Coulomb blockade regime is possible— ' 24 i 25 ' 26 and in 
agreement with the sequential-tunneling description. 



III. CURRENT AND FANO FACTOR 

The results of our calculation are summarized in Fig. [TJ 
which shows the average current (top panel) and the Fano 
factor (bottom panel) as a function of the bias voltage V 
and the gate voltage V g . Considering the current, the 
most visible feature is the Coulomb blockade diamond 
(green in Fig. [TJ. As long as the bias voltage is lower 
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FIG. 1: (Color online) Current (top) and Fano factor (bot- 
tom) in the one-electron blockade regime as a function of gate 
voltage V g and bias voltage V . The current is given in units 
of 67. Note the various plateaus in the Fano factor within 
the blockade regime, which are invisible in the current. The 
letters label the subregions as referred to in the main text. 
The dashed black line in the bottom plot indicates the line 
along which Figs. [3] and 2] are plotted. In this plot we use 
the parameters uh = 17to, W12 = 10to, 7l = Jr = 7 and 
k B T = O.Olto (and hence J = 0.24t and <f> = 0.79). 



than the energy difference between the energetically low- 
est state and the next state with one electron more or 
less, i.e., eV < Ei(V g ) — Ej(V g ) with i and j being two 
states differing by one electron, electron tunneling is ex- 
ponentially suppressed. In our model, we have three such 
blockade regimes with 0,1, and 2 electrons in the dot for 
zero bias voltage. The two parts of Fig. Q] show the cur- 
rent and the Fano factor for the most interesting of these 



blockade regions, the one-electron blockade regime. We 
will concentrate our discussion on this region, since most 
interesting features can be discussed here. 




-15P 

A-, — 1 C 1 B 11 — A, 

F D 1 



FIG. 2: (Color online) Energy of the states as a function of 
gate voltage V g in units of to. The capital letters correspond 
to the regions in Fig. [T] The parameters are chosen like in 

Fig.m 

We first look at the energies of the states in the block- 
ade regime as a function of the gate voltage, see Fig. [2] 
The capital letters in Figs. Q] and correspond to these 
blockade regions, and we will now show that the features 
in the Fano factor can be explained by the energy plot 
for the states in this region. We first consider the re- 
gions A. They are both characterized by a Fano factor 
F = 1, typical tunneling processes through a simple two- 
level system. Indeed, Fig. [5] shows that the first excited 
states differ in electron number. In A\ the excited state 
is |0) and in A 2 \S). Lowering the gate voltage, we come 
to region B, where the first excited state is the state |— ), 
i.e., a state with one electron as well. Thus, the descrip- 
tion in terms of a single state, which can be occupied 
or not, does not hold any more, and we observe an in- 
crease in the Fano factor to F = 2. In this region, the 
state |0) is still the first state available for tunneling, but 
the blocked state |— ) is energetically lower. Continuing 
to region C, in which the energies of the two-electron 
states become lower than the zero-electron state. So the 
first available excited state is the singlet The dif- 
ference between the regions C and D is somewhat more 
complicated. Since the difference Eg — E- is smaller 
than Es — E+, the state |+) still blocks the dot, if it 
is occupied. After a thermal excitation to the state \S), 
electrons tunnel through the states 15) and |— ) until an 
electron occupies the state |+) again and blocks the dot. 
In the regions C and D the Fano factor depends on the 
parameter <fi through the prefactor of the singlet states 
(see Eq. @). The ^-dependence is illustrated in Fig.[3]by 
varying the parameter un/to- As can be seen the posi- 
tions of the steps change simultaneously with the plateau 
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height. However, the height depends only on the interac- 
tion parameter, as we have checked by numerically vary- 
ing uh and </> independently. We will derive below ana- 
lytic expressions for the Fano factor, which confirm this 
behavior. 

The effect of temperature will be to wash out the sharp 
steps and plateaus. This is illustrated in Fig. |4] In- 
creasing the temperature to a value of the order of the 
level splitting first the sharp steps are washed out. The 
plateaus C and D, however, remains still visible and the 
values are still given, at least approximately, by the low 
temperature results. Finally, the temperature leads to a 
vanishing of the plateaus, which make it difficult to ex- 
tract the interaction parameter </> from the Fano factor. 
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FIG. 3: (Color online) Fano factor for different values of 
(j) in the regions C and D as a function of gate voltage 
(in units of to/e). The bias voltage is eV = 2to and 
ksT — O.Olto- Furthermore, the ratio U12/MH = 0.6 
is kept fixed, such that the different plots can be easily 
achieved in experiment by changing the interdot tunneling 
rate to. The parameters are chosen as (<f>\ J/to, 1112/to) — 
(0.62; 0.5, 5), (0.71; 0.34, 7), (0.79; 0.24, 10), (0.89; 0.12, 20) 
from right to left. The plateau heights, indicated by the 
black (region C) and gray (region D) dots, follow from the 
analytic expressions Eqs. {1} and J7J. 



IV. COUNTING STATISTICS 

A more detailed view of the transport characteristic 
can be obtained by looking at the full counting statis- 
tics (FCS)S1. The aim is here to find the cumulant gen- 
crating function (CGF) S(x), which is related to the 
probability P{N) of N charges passing through the sys- 
tem during the measurement time tM by exp5(x) = 
J2n P{N)e lNx . The knowledge of S(\) is equivalent to 
the knowledge of all the cumulants of the system accord- 
ing to Ck = (—id/d x ) k S(x)\ _ , such as current (first 
cumulant) and noise (second cumulant), I = eCxftu 



and S — 2e Ca/tM! respectively (where e is the elec- 
tron charge) . The simplest way to evaluate the counting 
statistics in this case is to use the method described in 
Ref. l2l|. We choose (without loss of generality) to count 
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FIG. 4: (Color online) Fano factor as a function of gate volt- 
age (in units of to/e) at various temperatures for <j> = 0.79, 
u h = 17io, and u\2 = 10io- The steps in the Fano factor are 
washed out by the temperature, however the noise remains 
super-Poissonian in the whole regime shown. 



the charges in the left lead. Therefore, we have to change 
the left tunneling rates in the off-diagonal elements of the 
transport matrix M: — > r^exp^x). The smallest 
eigenvalue of M, which we denote by Ao(x), determines 
the CGF as S( X ) = -t M Xo(x)- 

If we restrict our calculations to a certain region and 
thus reduce the number of involved states, it is possible 
to get an analytical expression for the Fano factor—. For 
region C, we include the one-electron states |+), |— ), and 
the two-electron singlet \S). Here, the relevant energy 
differences are both smaller than the bias: Eg — E- < 
Es — E + < eV/2. Therefore both Fermi functions for the 
tunneling processes are exponentially suppressed, but one 
is much bigger than the other: 

1 » f(E s -E+- eV/2) » f(E s — -EL — eV/2). (4) 

Therefore, we neglect f(Es — E— — eV/2) and take the 
rates to lowest order in the parameter x = f{E$ — E + — 
eV/2). Then the tunneling rates are 



r 



1 



+s 



1LX, 



l + < 



1 + 2 

2cf> 2 



1 



-(7£ + 7b)- 



Here, ^l/r are the bare tunneling rates for the left/right 
leads. Without loss of generality, we will count the 
charges in the left lead, which means, we have to replace 
7l 1l exp(±zx) (the sign depending on whether an 
electron enters or leaves the dot respectively) in the off- 
diagonal elements of the matrix M. The counting statis- 
tics is obtained from the smallest eigenvalue Ao, which 
we determine to lowest order in x. The result is 



s(x) 



7l7.r(1 + 4> 2 ) 



1 



7i?(l + 2 ) -1l 1 -p{4>)e>x 



(5) 
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where p(<f>) = 2 7r/[7l + 7r(1 + 2 )]- Finally we obtain 
the Fano factor 



F(<f>) 



7l + 7R 



(6) 



which is independent of the bias voltage as long as we 
are in the one-electron Coulomb blockade regime. For a 
symmetric structure with 7^ = 7^ we obtain Eq. ([1]). 

A similar result may be obtained for the region D in 
Fig. Q] Here, the situation is a bit more complicated, 
since the lowest-lying triplet state also plays a role and, 
thus, the problem involves four states. In the same way as 
before we neglect J(Et — E± — eV/2) since the transition 
between \S) <-> |+) is the least suppressed. Assuming 
that the temperature is much smaller than the singlet- 
triplet splitting, we calculate the series expansion for the 
Fano factor (with symmetric leads) to be 



F{4>) = 



2 + 170 2 + 32<?(> 4 + 160 6 



(7) 



2 + 110 2 + 16</> 4 + 3</> 6 
Both these results, Eqs. (JTJ) and 0, are in agreement 



with our numerical calculations as shown in Fig. [3j 



V. CONCLUSIONS 

In conclusion, we have shown that noise measurements 
in the Coulomb blockade regime of a double quantum 
dot can reveal interesting information about the quan- 
tum correlations created during tunneling. In particular, 
we have shown that the Fano factor in the one-electron 
blockade region is super-Poissonian and can be used to 
determine the interaction parameter (f>. This allows to 
measure the degree of entanglement of two electrons in 
the double quantum dot. 
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